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Voltage-controlled  switching  in  lateral  VO2  nano-gap  junctions  with  different  gap  lengths 
and  thermal  properties  was  investigated.  The  effect  of  Joule  heating  on  the  phase  transition 
was  found  to  be  strongly  influenced  by  the  device  geometry,  the  contact  material,  and  the 
current.  Our  results  indicate  that  the  V02  phase  transition  was  likely  initiated  electroni¬ 
cally,  which  was  sometimes  followed  by  a  secondary  thermally-induced  transition. 
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The  metal-insulator  phase  transition  in  the  correlated  electron  material1,  vanadium  dioxide 
(V02),  has  attracted  significant  attention  in  optics2-9  and  electronics10-15  for  its  applications  in  low 
power,  compact  and  high  speed  switching.  This  is  because  the  conductivity  of  the  V02  can  change 
by  up  to  5  orders  of  magnitude  in  response  to  external  stimuli  such  as  heating  beyond  340  K16, 
applied  electric  fields17,18,  surface  charge  accumulations11,  mechanical  strain19,  and  terahertz20  or 
optical  pulses21.  Additionally,  a  main  advantage  of  V02  is  that  a  phase  transition  based  on  elec¬ 
tronic  reconfiguration  can  occur  at  picosecond  time  scales21-23. 

Many  recent  experiments  have  investigated  the  electronic  phase  transition  in  lateral24-27  and 
vertical28,29  V02  junctions.  However,  most  time  domain  studies  have  shown  switching  times  that 
are  at  best  tens  of  nanoseconds28,30,31,  which  suggest  that  the  phase  transition  may  be  limited  by 
thermal  transients32.  Although  the  possibility  of  a  thermally  induced  phase  transition  is  ruled  out 
by  heat  transfer  simulations  of  single  crystalline31-33  and  polycrystalline  V02  junctions34,  exper¬ 
iments  indirectly  measuring  the  local  temperature  of  two  terminal  V02  junctions  suggest  that  the 
temperature  rose  above  the  insulator-metal  phase  transition  temperature  of  V02  at  Tj^m  =  340 
K35,36.  These  observations  suggest  that  the  role  of  thermal  effects  on  the  phase  transition  can  vary 
significantly  based  on  the  parameters  of  the  experiment.  As  a  result,  understanding  these  param¬ 
eters  not  only  sheds  light  on  the  nature  of  the  phase  transition,  but  may  also  help  improve  the 
electronically-induced  switching  times  by  showing  how  thermal  effects  may  be  suppressed. 

In  this  work,  we  study  the  voltage-controlled  switching  of  V02  lateral  junctions  with  different 
lengths  and  contact  metals  to  determine  the  electric  field  at  the  onset  of  phase  transition  at  various 
heat  dissipation  rates.  We  show  that  when  the  current  is  low,  the  electric  field  at  the  onset  of  the 
insulator-metal  transition  is  independent  of  the  thermal  conductivity  and  work  function  of  the  con¬ 
tact  metal,  but  the  reverse  metal-insulator  transition  and  hysteresis  width  are  strongly  dependent 
on  the  contact  metal.  Our  results  suggest  that  the  phase  transition  is  initiated  electronically,  which 
is  sometimes  accompanied  by  secondary  thermal  effects.  The  thermal  contribution  is  determined 
by  the  device  geometry,  the  contact  metals,  and  the  external  load  resistance.  The  variations  of 
these  parameters  in  different  reports26,32,35'36  help  explain  the  discrepancy  in  the  observed  cause 
of  the  phase  transition.  Our  results  indicate  that  thermal  effects  may  be  reduced  or  potentially 
suppressed  through  device  design. 

We  fabricated  a  series  of  lateral  V02  junctions  as  illustrated  in  Fig.  1(a).  At  —  100  nm  V02 
film  was  deposited  using  magnetron  sputtering  of  a  vanadium  target  on  a  2  //m  thick  thermally 
grown  Si02  on  a  Si  substrate3.  V02  lateral  junctions  were  fonned  by  electron  beam  lithography 
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FIG.  1.  (a)  The  top  and  side  view  schematic  diagram  of  a  VO2  lateral  junction,  (b)  A  SEM  of  the  device 
and  the  magnified  image  of  the  gap.  (c)  The  simulated  electric  field  profile  and  (d)  current  density  profile 
when  3  V  was  applied  over  a  300  nm  gap. 


and  lift-off  of  two  metal  contacts  with  a  width  of  W  =  10  /mi,  thickness  of  100  nm,  and  a 
separation  gap  of  length  L.  Figure  1(b)  shows  a  scanning  electron  micrograph  (SEM)  of  a  device 
with  a  gap  of  300  nm  and  palladium  (Pd)  contacts.  The  inset  shows  a  magnified  image  of  the  gap. 

When  a  voltage  is  applied  between  the  two  electrodes,  an  electric  field  establishes  across  the 
gap  and  a  small  current  leaks  through  the  VO2  film.  Figure  1(c)  shows  the  simulated  electric  field 
profile  when  3  V  is  applied  across  a  300  nm  gap.  Although  a  large  portion  of  the  field  resides 
above  the  film  in  the  air  gap,  the  average  field  inside  the  V02  is  75  kV/cm.  The  field  strength 
is  consistent  with  the  reported  field  values  required  to  induce  an  electronic  phase  transition  in 
V0228’30.  Figure  1(d)  shows  the  corresponding  current  density  profile  of  Fig.  1(c)  inside  the  V02 
slab,  which  is  mainly  concentrated  between  the  two  electrodes. 

Figure  2(a)  shows  a  simplified  circuit  model  of  the  device  and  experiment.  The  VO2  junction 
(shaded  area)  is  modeled  as  a  capacitance,  C(,ap,  in  parallel  with  a  resistance,  Fl(jap.  A  voltage,  V'app, 
was  applied  across  each  device  and  an  external  resistor  RL  that  limited  the  current.  The  current,  I, 
and  the  voltage  drop  across  the  gap,  Vr(iap,  were  monitored.  At  a  sufficiently  high  applied  voltage, 
the  insulator-metal  transition  would  occur  in  the  V02,  resulting  in  a  drop  in  the  resistivity  and  an 
increase  in  the  current. 

Figure  2(b)  shows  /  vs.  Vdap  as  Vapp  was  ramped  up  (green  line)  and  down  (purple  line)  for  a 
device  with  Pd  contacts,  L  =  300  nm  and  /fi  =  550  Q.  The  numbered  arrows  in  the  figure  trace 
Voap-  For  small  values  of  Vapp  at  point  (1),  the  V02  was  in  the  high  resistivity  phase,  so  the  current 
increased  slowly  with  voltage.  At  point  (2),  the  voltage  increased  beyond  the  insulator-metal  phase 
transition  voltage,  so  the  resistivity  of  the  V02  reduced,  resulting  in  a  sharp  increase  in  the 
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FIG.  2.  (a)  The  simplified  circuit  model  of  the  device  and  experiment.  The  (b)  voltage-current, (c)  voltage- 
resistivity,  and  (d)  voltage-power  relationship  of  a  junction  with  a  gap  length  of  300  nm  and  Pd  contacts. 

current  accompanied  by  a  drop  in  Vbap,  or  a  large  negative  differential  resistance.  After  point  (3), 
the  VO2  was  in  a  low  resistivity  phase,  and  the  current  increased  linearly  with  the  applied  voltage. 
When  V'app  was  reduced  back  to  zero  between  points  (5)  to  (7),  the  curve  traced  a  different  path 
due  to  the  hysteresis  of  V02.  As  VApp  was  reduced,  VGap  increased  until  point  (6)  when  it  reached 
the  metal-insulator  transition  voltage,  Vm^i,  after  which  the  V02  returned  to  its  high  resistivity 
phase. 

Figure  2(c)  shows  the  resistivity  of  the  underlying  V02  film  calculated  as  p  —  RGapWt/L 
when  Rl  =  550  Q  and  10  kQ.  For  both  values  of  RL,  p  was  relatively  constant  for  small  values 
of  Vdap,  but  it  dropped  significantly  when  Vbap  exceeded  and  increased  to  its  original  value 

when  Vbap  was  reduced  below  Vm^i-  As  expected,  the  initial  value  of  p  was  independent  of  RL, 
but  Vi-yM  and  the  magnitude  of  the  drop  in  p  were  dependent  on  RL.  As  we  shall  show,  this 
dependence  stems  from  different  Joule  heating  of  the  V02  slab37. 

To  examine  the  Joule  heating  effects  in  our  devices,  we  calculated  the  dissipated  electrical 
power  in  the  device  by  P  =  /  x  VGap.  Taking  P  to  be  the  Joule  heat  source,  we  simulated  the 
expected  temperature  in  V02.  Figure  2(d)  shows  P  (left  axis)  and  the  simulated  temperature36,  T, 
(right  axis)  of  the  V02  slab  as  a  function  of  Vbap-  The  simulated  temperature  is  below  Tf^M  =  340 
K  at  the  onset  of  the  phase  transition  for  both  values  of  f?L,  indicating  that  Joule  heating  alone 
would  not  be  sufficient  to  initiate  the  phase  transition.  However,  for  RL  =  550  Q  and  at  any 
given  Vbap,  the  simulated  temperature  in  V02  is  higher  than  that  of  R\_  =  10  1<Q  which  would 
lower  Vj^m-  This  temperature-induced  decrease  in  Vj^M  has  been  observed  in  bulk  V02  films 
and  has  been  studied  by  many  groups  for  the  electrical27’34,36’38,  optical39  and  strain  induced  phase 
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transitions19.  For  our  devices,  we  found  that  when  RL  was  10  kQ,  an  increase  in  temperature  of  8 
K  would  result  in  the  same  VI^fM  as  the  case  when  Rh  =  550  Q.  This  value  agrees  well  with  the 
simulated  temperature  increase  of  6.5  K,  as  shown  in  Fig.  2(d). 

Joule  heating  can  also  explain  the  difference  in  the  magnitudes  of  the  resistivity  drop  after 
the  phase  transition  for  the  two  values  of  J?L.  P  at  the  onset  of  the  transition  was  similar  for 
both  values  of  /?L,  indicating  that  the  same  energy  was  required  to  initiate  the  phase  transition. 
However,  the  change  in  P,  and  consequently  T,  after  the  transition  was  negative  when  RL  =  10 
kO  and  positive  when  R\  =  550  Q.  For  the  case  of  ll\_  =  550  O,  the  increase  in  P  after  the 
transition  would  lead  to  a  rise  in  temperature  of  the  V02  slab  above  T/_>m  =  340  K,  which  could 
initiate  a  secondary  thermally-induced  phase  transition.  The  thermal  transition  has  been  shown  to 
exhibit  a  larger  resistivity  switching  magnitude  than  the  electronic  transitions28.  Thus,  a  thermal 
transition  would  account  for  the  larger  change  in  p.  This  change  in  p  was  consistent  with  thermal 
resistivity  switching  of  our  bulk  films.  When  RL  =  10  kQ,  however,  the  decrease  in  P  after  the 
transition  would  lead  to  a  decrease  in  the  temperature,  which  would  account  for  a  smaller  change 
in  p. 

The  simulated  temperatures  before  and  after  the  phase  transition  for  the  two  values  of  RL 
suggest  that  the  phase  transition  was  most  likely  initiated  electronically  and  was  accompanied  by 
secondary  thermal  effects.  To  verify  that  the  phase  transition  was  not  initiated  by  Joule  heating,  we 
fabricated  four  sets  of  devices  with  different  thennal  dissipation  rates  using  metals  with  different 
thermal  conductivities  as  the  contacts.  This  allowed  us  to  vary  the  heat  dissipation  rate  of  the 
devices  while  keeping  the  heat  generation  rate  constant.  Since  different  contact  metals  would 
lead  to  different  work  functions  mismatches  and  V02-metal  interface  resistances,  we  swept  L  to 
extrapolate  the  electric  field  at  the  onset  of  the  transition  and  remove  any  work  function  differences 
at  the  metal-V02  interfaces. 

DC  measurements  similar  to  those  of  Fig.  2  were  conducted  on  samples  with  Pd,  gold  (Au), 
silver  (Ag),  or  copper  (Cu)  contacts.  All  samples  were  made  from  the  same  V02  film.  They  all 
had  a  width  of  W  =  10  //m,  a  metal  thickness  of  100  nm  and  a  separation  gap  of  length  L,  which 
varied  between  150  nm  to  1000  nm.  RL  was  fixed  at  10  kQ.  Figure  3(a)  shows  the  measured  Vj^m 
for  various  values  of  L  (scattered  points)  and  the  corresponding  linear  fits  (dashed  lines)  for  each 
contact  metal.  The  error  bars  represent  the  standard  deviation  in  Vi^m  for  multiple  measurements. 
The  y-intercepts  of  Fig.  3(a)  (V/_^m  as  L  — >•  0)  are  related  to  the  work  function  difference  at  the 
metal-V02  interface  due  to  the  interface  resistance  and  the  work  function  difference  between  the 
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V02  and  the  contact  metal.  The  intercepts  of  these  plots  are  2.2  ±  0.2  V  for  Pd  ,  2.3  ±  0.3  V  for 
Au,  3.3  ±  0.5  V  for  Ag,  and  3.2  ±  0.3  V  for  Cu.  The  order  of  these  values  correlates  with  the 
work  function  differences  between  the  metal  and  that  of  VO2  reported  in  Ref.40.  After  the  phase 
transition,  these  interface  effects  decrease  due  to  both  the  change  in  the  V02  work  function  and 
the  increase  in  carrier  density.  These  measurements  also  show  that  neglecting  interface  effects 
results  in  an  overestimation  of  the  transition  electric  field  by  up  to  63  %  for  a  gap  length  of  300 
mn.  Although  these  interface  effects  may  be  negligible  for  L  >  1  //m,  they  can  lead  to  significant 
power  dissipation  in  nanometer  scale  junctions.  Doping  may  be  used  to  create  ohmic  contacts  to 
suppress  these  interface  effects. 


L(nm) 


Thermal  Conductivity  (W/m.K) 


FIG.  3.  (a)  The  forward  transition  voltage  ( )  of  devices  with  different  thermal  dissipation  rates,  (b) 
The  electric  field  at  the  onset  of  phase  transition  for  the  forward  and  reverse  transitions  as  the  function  of 
metal  thermal  conductivity. 


As  a  measure  of  the  “average”  electric  fields  in  the  gap  for  the  forward  (E^m)  and  reverse 
(Em^j)  phase  transitions,  we  calculate  the  slopes  of  V^M  and  VM~>i  vs.  L.  Figure  3(b)  shows 
E^m  and  EM_,[  as  a  function  of  the  thermal  conductivity  of  the  contact  metal41.  We  found 
that  the  values  of  E^M  were  similar  for  the  different  contact  metals  despite  their  vastly  different 
thermal  conductivities,  indicating  that  the  forward  transition  was  unlikely  to  be  thermally-induced. 
If  the  transition  was  thermally  initiated,  the  Pd  devices  with  the  lowest  thennal  conductivity  would 
have  required  the  smallest  value  of  Ej^m  for  the  phase  transition,  which  was  clearly  not  the  case. 

On  the  other  hand,  Em^j  was  strongly  dependent  on  the  metal  thennal  conductivity,  which  can 
be  explained  by  the  dissipated  power  and  temperature  of  the  V02  after  the  phase  transition  as  in 
Fig.  2(d).  For  RL  =  10  kfl,  the  dissipated  power  decreased  at  the  onset  of  the  phase  transition. 
Well  after  the  phase  transition,  however,  the  metallic  V02  acted  as  a  resistive  heating  element 
that  raised  the  temperature  with  increasing  voltage.  At  a  sufficiently  high  power  density,  the  local 
temperature  of  V02  could  be  higher  than  the  metal-insulator  phase  transition  temperature,  Tm^-i 
=  317  K.  Consequently,  the  reverse  phase  transition  was  dominated  by  thermal  dissipation  rates, 
which  occurred  at  lower  EM^T  values  for  metals  with  lower  thermal  conductivities.  Therefore, 


6 


devices  that  used  contact  metals  with  higher  thermal  conductivities  exhibited  smaller  electric  field 
hysteresis  widths.  It  should  be  noted  that  the  value  of  Em^i  is  generally  lower  than  Er_tM,  due  to 
the  hysteresis  of  VO2  and  the  change  in  thermal  conductivity  of  VO2  after  the  phase  transition.  The 
temperature  dependence  of  the  electronic  transition  explains  the  observed  dependence  of  Vj^M  on 
the  external  load  as  shown  in  this  work  and  similarly  in  Ref.37.  For  small  values  of  /fi,  compared 
to  the  V02  channel  resistance,  the  jump  in  the  power  density  at  the  onset  of  the  phase  transition  can 
substantially  raise  the  temperature,  resulting  in  strong  thennal  contributions  to  the  phase  transition. 


TABLE  I.  Comparison  of  thermal  effects  in  the  voltage  switching  of  recent  two  terminal  VO2  switches. 


No. 

VO2  Film  Properties 

Device  Geometry 

Properties  at  I  — >•  M  Phase  transition 

Deposition 

Substrate 

Film 

p  (fl.cm) 

p  ratio 

L 

W 

R  L 

Contact 

{E^mT 

Power 

Temperature 

Technique 

thickness 

(insulating 

Max(p) 

{pm) 

t/'m) 

(kO) 

(kV/cm) 

density 

increase  (K) 

(nm) 

phase) 

Min(p) 

(pW/pm3) 

Pd 

21° 

0.42 

17 

Sputtering 

Si02 

10 

Au 

23c 

0.46 

12 

lb 

100 

2 

100 

0.5 

10 

Cu 

23c 

0.46 

11 

0.5 

Pd 

20c 

0.6 

18 

2d 

Sputtering 

A1 

370 

- 

146 

0.37 

300 

- 

Ti/Au 

50 

5.5xl0-5 

- 

3e 

Sputtering 

Si02 

200 

0.2 

80 

16 

8 

0.2 

Au 

5 

1.875 

77/ 

4 9 

Sputtering 

A1203 

130 

4 

1000 

20 

50 

180 

Au 

20 

0.2 

17  f 

5h 

Sole  Gel 

ai2o3 

100 

0.2 

200 

10 

10 

15 

Ni 

10 

0.5 

27 

6l 

Pulsed  Laser 

ai2o3 

100 

- 

- 

25 

30J 

1 

Pt 

12 

4 

180 

Deposition 

lk 

Evaporation 

ai2o3 

90 

0.002 

400 

0.15 

10 

0.05 

Ti/Au 

9 

0.5 

- 

of  V 

a  The  Vi^m/L  value,  unless  otherwise  specified.  6  This  work.  c  Extrapolated  from  the  slope  of  Vi^m  vs.  L. 

d  Reference25.  e  Reference33.  f  Experimentally  measured  value. 

9  Reference32.  h  Reference30.  1  Reference34. 

J  Width  of  the  observed  filament.  k  Reference27. 


Table  1  summarizes  the  VO2  film  properties,  the  device  geometry  and  the  estimated  contri¬ 
bution  of  electrical  and  thermal  effects  in  this  work  and  other  recent  experiments  on  electrical 
switching  of  V02.  The  properties  of  V02  are  also  provided  for  comparison,  where  f>\  is  the  re¬ 
sistivity  of  the  insulating  phase  and  p  ratio  is  Max(p)/Min(p).  The  device  geometry  and  external 
load  resistance  have  strong  effects  on  the  temperature  increase  in  these  devices.  The  reported 
temperature  increase  is  directly  proportional  to  the  power  density  and  provides  a  metric  for  com¬ 
parison.  Devices  with  small  cross-sections  and  low  loads  have  high  dissipated  power  densities 
and  consequently  large  temperature  increases.  Similarly,  in  our  work,  when  RL  was  reduced,  the 
temperature  at  the  onset  of  the  phase  transition  increased.  Devices  with  higher  temperatures  at  the 
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phase  transition  generally  required  a  lower  Er^M.  This  is  expected  as  increasing  the  temperature 
increases  the  carrier  density  inside  the  V02  and  consequently  reduces  Em^j3H.  It  should  be  noted 
that  devices  (3)  and  (6)  had  small  cross-section  areas  and  long  gaps,  which  created  sufficiently 
high  power  densities  to  increase  the  V02  temperature  beyond  Tj^M-  In  these  cases,  the  observed 
phase  transitions  were  mostly  thermally-induced. 

In  our  experiments,  although  the  transition  was  initiated  electrically,  it  remains  unclear  whether 
the  phenomenon  was  a  field17,  a  current  or  an  avalanche  effect34.  Our  preliminary  experiments  on 
electric  field  switching  in  similar  structures,  where  the  current  was  completely  suppressed  using  a 
thin  high-k  dielectric  that  separated  the  contacts  and  the  V02,  showed  that  a  purely  electric-field 
induced  switching  could  not  be  achieved  at  similar  field  strengths.  The  high-k  dielectric  typically 
experienced  catastrophic  breakdown  before  the  V02  phase  transition  occurred. 

In  conclusion,  we  have  demonstrated  that  voltage  controlled  switching  in  two  tenninal  V02 
switches  in  the  nano-scale  is  dominated  by  electronic  effects  if  the  current  could  be  suppressed. 
A  large  current  causes  a  rapid  increase  in  the  power  density  after  the  phase  transition  to  result 
in  a  secondary  thermal  transition.  High  frequency  circuits  for  current  regulation  after  the  phase 
transition  may  be  a  solution  to  suppress  secondary  thermal  effects.  These  circuits  are  typically 
used  for  current  quenching  in  avalanche  photo-diodes42  and  can  be  adapted  to  V02  devices  to 
achieve  picosecond  switching  times22,23. 
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